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Abstract 



In this paper we analyze the process in a model 11 2HDM and MSSM. All the leading 

terms of Wilson coefficients relevant to the process are given in the large tan/3 limit. It is shown that 
the decay width for Bs depends on all parameters except m^o in the 2HDM. The branching 

^ ratio of Bs — >■ can reach its experimental bound in some large tan/3 regions of the parameter 

O . space in MSSM because the amplitude increases like tan^/? in the regions. For l=r, the branching 
^ [ ratio can even reach 10~^ in the regions. Therefore, the experimental measurements of leptonic 
\^ decays of B^ could put a constraint on the contributions of neutral Higgs bosons and consequently 
O ■ the parameter space in MSSM. 
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1 Introduction 



Bg —* l^l^, as one of flavor changing neutral current processes, is sensitive to structure of the 
standard model (SM) and new physics beyond SM, and is expected to shed hght on the existence of 
new physics before the possible new particles are produced at colliders. Theoretically, the process 
is clean because only the nonperturbative quantity involved is of the decay constant of Bg and it 
is relatively easy to be calculated by so far well-known nonperturbative methods such as QCD sum 
rules, lattice gauge theory, Bethe-Salpeter approach, etc. Therefore, it provides a good window 
to probe new physics. Experimentally, the 95% confidence level upper bound on the Bg — > 
branching fraction has been given [jl|: 

Br{Bs fi+fi-) < 2.6 X 10"^ (1.1) 

The planned experiments at B- factories are likely to measure branching fractions as low as 1O^^[0]. 

Compared to the rare decay B Xg'-f, Bg —>■ (as well as B ^ XJ^l~) is of more advantage 
for the study of the Higgs sector in the large tan/5 case in a model II two Higgs doublet model 
(2HDM) or supersymmetric models (SUSY) since the contributions to -B — Xs'y coming from Higgs 
sector are indeed independent of tan/3 when tan/3 is larger than a few ( say, 4 ). 

The branching ratio for Bg has been calculated in SM and beyond SM in a number of 

papers p, ^, |[ . In a recent paper the process in a model II 2HDM with large tan/? is reanalyzed. 
It is correctly pointed out that the contributions of the box diagram to this decay at the leading 
order of tan/? are missed and a minus for the contribution of penguin diagram involving and 
in the loop is also missed in the earlier literature 0, ^, However, there are some points in 
the paper which need to be clarified. First, the argument that the trilinear H^H^H{ }l=h^,H^) 
couplings should not be considered as tan/? enhanced is not correct, as we shall argue below. Second, 
although the contribution of box diagram is the same order as those of penguin diagrams in the 
large tan/3 limit it is numerically smaller than those of penguin diagrams and consequently the claim 
that the box diagram gives the dominant contribution in the 't Hooft-Feynman gauge is not true. 
In the paper we shall give a detailed argument (as we know, there is no such argument presented in 
the literature ) to show why their claim on H^H^H couplings is not correct (see section 3). These 
arguments are important to clarify where are the disagreements in the literature and make one have 
a correct conclusion on the decay in 2HDM. 

There are more box diagrams in SUSY than that in 2HDM. The contributions of box diagrams in 
the analysis in supersymmetric models are missed in the previous papers [§, 0, § . The contributions 
are also omitted in the refs. |]14|, |16| since they neglect the mass of a lepton in calculating Wilson coef- 



ficients. However, for l=/i, r the contributions in the large tan/3 case are important and consequently 
should not be neglected. In the paper we calculate the contributions to Wilson coefficients Cg, Cg. 
from the box diagrams and carry out a complete analysis in a model II 2HDM and SUSY with large 
tan/?. 

The paper is organized as follows. Section 2 is devoted to the effective Hamiltonian responsible 
for h — > sl^l~ . We calculate the contributions to Wilson coefficients from the box diagrams and 
give all the leading terms of Wilson coefficients in a model II 2HDM and SUSY with large tan/3 in 
section 3. In section 4 we present the numerical results. In section 5 conclusions are drawn. Finally 
the contributions to Wilson coefficients Cg, Cq. from individual diagrams in a 2HDM and MSSM are 
given in the appendix. 
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2 Effective Hamiltonian for b si^i 

The effective Hamiltonian describing the ffavor changing processes b — > si'^i" can be defined as 

/-f 10 10 

H^jf = -^A,(EC.(/i)0,(/.) + ECQ.(/i)g.(/x)). (2.2) 

1=1 i=i 

where = VtbV^*, OjS (z = 1, ■ ■ ■ , 10) are the same as those given in the ref. |TD|, |T2|]Q, and QiS come 
from exchanging neutral Higgs bosons and have been given in refs. p, ^. 

The QCD corrections to coefficients Cj and Cq- can be incorporated in the standard way by 
using the renormalization group equations. Qi{i = 1, ■ ■ ■ , 10) does not mix with Og, Og so that the 



evolution of Cs and Cg remains unchanged and are given in ref.||10 



Att A ,8 

Cs{m,) = Csimw) + ^^^^h^i^ - r''"') + ^(1 - r'^'^'W^imw), 

= C^{mw) + 1.937C2(mvK) (2.3) 

Cg{m,) = Cg{mw). (2.4) 

It is obvious that operators Oiii = 1, ■ ■ ■ , 10) and Qi{i = 3, ■ ■ ■ , 10) do not mix into Qi and Q2 
and also there is no mixing between Qi and Q2. Therefore, the evolution of Cq-^,Cq^ is controlled 
by the anomalous dimensions of Qi,Q2 respectively. 

CQ^im,)=v-''^/^'CQXmw), ^ = 1,2, 

= 1.2ACQXmw) (2.5) 

where 7q = —4 Jill loop anomalous dimension of s^feij, r] = ^"l^j^^^ ~ 1.72, and /3o = 

11 - (2/3)n/ = 23/3. 

For the decay Bg —>■ , the matrix element of Tie// is to be taken between vacuum and \B^) 
state. Because 

(0| s a^"" Pr b = (2.6) 

and the Og term in eg. ( ^.2[ ) gives zero on contraction with the lepton bilinear due to = p+ + pi, 
only the operators Og and Qi 

Og = SLYbiHl^l, (2.7) 

Oq, = SLbJl, (2.8) 

Oq, = SLbnh^l (2.9) 

are involved and the important thing we need to do is to calculate the Wilson coefficients of the 
operators at /i = mw- Cg{my/) has been calculated in SM in a 2HDM [10| and in SUSY mod- 
els |I4| , respectively. Cg{mw) in the 2HDM is the same as that in SM for the large tan(3 scenario. 

^We follow the convention in ref. for the indices of operators as well as Wilson coefficients. In the convention 
of ref.|l|] Og (O9) is changed into O9 (Oio). 
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The box diagram contributions to Cg^niw) (as well as Csi^m^y)) in SUSY which are proportional 
to tan^/? are missed in the previous calculations |]16[ ||, H- Cq^s have also been calculated in 
2HDM [|, ^ and in SUSY 0, 0. However, some leading terms in the large tan/3 limit are missed 
in the previous papers. We shall calculate Cg and Cq. at /i = mw in the next section in order to 
give a complete and correct result. 

By using the equations of motion for quark fields, we have 



2 

m, 



{0 \s 1^ b\B',) = tfs^^^ (2.10) 

where /^^ is the decay constant of Bg defined by 

{0\sri5b\B',) = -tfB^p>^ (2.11) 
Thus the effective Hamiltonian ( p.2|) results in the following decay amplitude for Bs — > 

A = ^0§^'^BJB^t [ Cq.uv + {Cq, + 2miC,) u^,v] , (2.12) 

where = mi/mB^ and itlb, ^ mb + nis have been used. Then it is straightforward to obtain the 
branching ratio for Bs l^l^ 

^ F^EM 



Br (Bs - i^n = ^^^<r^J4|At| VI - ^mj 



(l - 4mj) ICgf + \Cq, + 2m,C9p] , (2.13) 



where tb is the B^ lifetime. 



3 Wilson coefficients 
3.1 In 2HDM 

Consider two complex hypercharge F = 1, SU{2)w doublet scalar fields, 0i and 02- The Higgs 
potential which spontaneously breaks SU{2) x U{1) down to U{1)em and conserves CP symmetry 
can be written in the following form [ p!7[| : 

V{<Pi, 02) = Ai(<l>I$i - vff + A2($^$2 - 

+A3[($l$i-t;?) + ($^2<^'2-t^2')F 

+ A4[($l$l)($^2'^'2) - (0l$2)($^$l)] 

+A5[Re($I<l>2)-t;it;2]' 

+A6[Im($l$2)]' (3.14) 
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Hermiticity requires that all parameters are real. If Aj > the potential is semi-positive and the 
minimum of the potential is at 



< $1 >= 



< $2 > = 



V2 



(3.15) 



thus breaking SU(2) x U{1) down to U{1)em- From the potential the mass eigenstates are easily 
found as follows. The charged Higgs states are 



G± = (j)f cos p + (f)2 sin p 
= —(pf sin P + (pf cos P, 
where the mixing angle /? is defined by tan/?=f2/fi. The CP-odd states are 

G° = ^2 [Jm0? cos p + Im(f)l sin p] 
A° = V2[-Im(pl sin p + Imcpl cos P] , 



(3.16) 



(3.17) 



The would-be Goldstone bosons and are eaten by the W and Z bosons. The physical Higgs 
boson masses are 



(3.18) 
(3.19) 



where f ^ = f ^ + f| is fixed by the W boson mass, = \g'^v'^. Diagonalizing the mass matrix of 
CP-even Higgs 



2 r Acos^P{Xi + A3) + sin'^pX^ (4A3 + X^)sin2p/2 
' [ (4A3 + A5)sm2/?/2 4sm2/5(A2 + A3) + cos^/^As 



(3.20) 



results in the CP-even eigenstates 



with masses 



— V2) sin a], 
= ^[{-Re(t)l - vi) sin a + {Recpl - V2) cos a] 



(3.21) 

(3.22) 
(3.23) 

A_cos2/3 - 4A2sm2/? + AXicos'^P ^^'^^^ 

where X± = 4A3 ± A5. In the potential ( |3.14D there are 8 parameters : Aj, i=l,...,6, vi and V2- 
Because, as said above, v"^ is fixed by mw, there are seven independent parameters in a general CP 



2 I 2 
rUf^o + m^o 



and mixing angle 



tan2a 



v'^[4:{cos'^pXi + sin^pX2) + A+] , 
v^[X_cos2P - AX2sin^P + AXiCos^ P]/ cos2a 

X+sin2P 
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invariant 2HDM. Six of them can be expressed in terms of mixing angles a and (5 and Higgs masses 
mH±,mjio,mjjo, 171)^0. The seventh needs to be fixed by measuring one of the quartic couphng in 
( p.l4| ). For simphcity we shall assume Ai = A2 hereafter so that we have six independent parameters 



in the model. Taking Ai = A2, eg s . ( |3 . 2 2| , |3 . 2 3| , |3l2^) reduce to 

mjjo + mlo = v^{AXi + A+) , (3.25) 



mjjo - mlo = vy{AXi + X^ycos^2(3 + Xlsin'^2(3 (3.26) 

tan2a = \ tan2(3. (3.27) 
A_ + 4Ai 

The three equations show explicitly that the angle a as well as masses mHo^m^o can be traded 
for Aj, i=l,3,5 (or, equivalently, 1,+,-) no matter how large tan/? is. That is, a, as one of the set 
of six independent parameters which contains both a and j3 as well as others, can take any value 
independent of tan/?, as it should be. Therefore, the statement in ref. that the angle a depends 
on (3 is not correct. When tan/3 approaches to infinity, if (4Ai + A_) and consequently m|^o — mlo is 
of order cot/?, say, (4Ai + A_)= c tan2/? with c a constant of order one., then sin2a is of order one. 
If (4Ai + A_) is of order one, then in the large tan/? limit it follows that sina ~ cotf3 or l-cot^fi /2 
so that sin2Q; goes always as cot/5 which cancels the tan/? enhancement. However, the conclusion is 
valid only at tree level. Once the radiative corrections are included it would change, which is similar 
to the situation that happens in the Higgs sector of MSSM, i.e., the radiative corrections violate the 
tree level mass relations and one treats Higgs boson masses as free parameters to be determined by 
experiments. Therefore, even in this case we should still treat a as well as (3 and Higgs masses as free 
parameters in the general 2HDM defined above so that the tan^/? enhancement due to the trilinear 
H^H^H{ H=/i°,if°) couplings should be considered, as we did in ref. ||^. 

As usual, in the model II 2HDM the Higgs-fermion Yukawa couplings are given by 



^Yuk 



-YdQ<^iD -YuQ^U -YiW + h.c. (3.28) 



where = 2x2$*. So down-type quarks and charged leptons (up-type quarks) acquire masses by 
their couplings to $1 ($2)- 

Feynman rules in the above general 2HDM have been given. Vertices with one or two gauge 



bosons and vertices involving two fermions and one boson are given in ref. |T8[. The three Higgs 
boson vertices can be found in ref. @]. The vertices involving one Goldstone boson and two Higgs 
bosons have also been given in ref. [113. We use these Feynman rules in calculations of Wilson 



coefficients. 

As pointed out in section II, for large tanjS, Cg{mw) in the 2HDM is the same as that in SM. 
The leading contributions to Cq- in the large tan/3 limit come from the diagrams in Fig. 1. In our 
previous paper we paid attention to the contributions of neutral Higgs bosons and missed the 
contribution from the box diagram involving one charged Higgs and one W boson which is of order 
tan^/3 in the large tan/3 limit |0 . We carry out a calculation for the diagram and confirm the result 
in ref. 0. In this paper we include the contribution and correct a sign for penguin. In order to 
separate contributions from individual diagrams, we write Cq- as 

CQ^ = C'+Cl+Cg, (3.29) 
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where Cq_, Cq., and Cq_ denote the contributions from self-energy type diagrams, box diagrams, 
and Higgs penguin diagrams, respectively. In appendix A we present all contributions proportional 
to tan^/3 in Feynman-t'Hooft gauge. Adding all tan^/5 contributions together, we have 



where 



facVt 



Inyt sin^{2a) {m\a - m^jof 



' facUt 



2 2 



fliVt 



i-yt 



nibe 



nibe 
rrirg 

0, i 



5, ■■■,10 



fa 



m 



yt 



m 



fiiv) 



l-y + ylny 



(3.30) 

(3.31) 

(3.32) 

(3.33) 
(3.34) 



The difference between eq.( |3.30| ) and the result in ref. is that the first term in the brackets in 
eq. (|3.30| ) is incorrectly omitted in ref. 0. It is worth to note that in the above equations mb=mi,{mw) ■ 



3.2 In SUSY 

In the minimal supersymmetric standard model (MSSM) or supergravity model (SUGRA) the Higgs 
sector is the same as that in a model II 2HDM by imposing the following constraints on the param- 



eters |ir8[| : 



A2 = Ai (3.35) 

\, = \{g^ + g^)-\, (3.36) 

A4 = 2Ai - ^^'^ (3.37) 

A5 = A6 = -2(Ai + 2A3). (3.38) 



And all Feynman rules can be found in ref. In addition to Fig. 1, the diagrams in Fig. 2 also 
give the leading contributions. Besides box diagrams, five different sets of contributions to the decay 
h — >■ sl^l~ are present in SUSY. They can be classified according to the virtual particles exchanged 
in the loop: a) the SM contribution with exchange of W~ and up-quarks; h) the charged Higgs boson 
contribution with H~ and up-quarks; c) the chargino contribution with x~ and up-squarks (m); d) 
the gluino contribution with g and down-squarks (c?); and finally e) the neutralino contribution with 
X° and down-squarks. As pointed out in refs. ]T^ , |T3|, |^, contributions from neutrilino-down type 
squark (e) and gluino-down type squark (d) loop diagrams are negligible compared to those from 
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chargino-up type squark diagrams because the flavor mixings between the third and the other two 
generations are small in minimal supergravity and constrained MSSM. Therefore, in addition to the 
SM (a) and charged Higgs (b) contributions, we only include the contributions from chargino-up 
type squark (c) loop diagrams in the paper. 

In some regions of the parameter space the dominant contribution to Cq. is proportional to 
tan^/3 and comes from the self-energy type diagrams, as pointed out in ref. [|], |^, |19|. In quite a large 
region of the parameter space the dominant contribution is proportional to tan^/3. Box diagrams can 
contribute terms with tan^/3 to Cq. which are greatly suppressed by {rUs/Mw)'^, therefore the largest 
contributions to Cq. from SUSY box diagrams remain proportional to tan^/?. Among the diagrams 
in Fig. 2 only the box diagram with charginos in the loop can give the tan^/3 enhancement to Cg{m\y). 
The contributions of the self energy type and penguin diagrams to Cq. have been calculated by us 
in refs. p, |^. We calculate the contributions of the box diagrams and summarize all contributions 
in the appendix B. Adding all contributions given in the appendix B, one has 

2 2m^om^o 



fco{xH-,l, 0)} - fac Y' 51 ^ibk'Ask\^ij^kk'rhH.lx~fBO [x - 



\f2 



ks ^ 



rhH^Jx~fBO (x^-, xsfe) + 0{tan^f3), (3.40) 

Cq^{Mw) = fac {Xt - l)fco{xH-, l,Xt) + /cO (xj/- , 1, 0) 

M, 2 6 



^b'^t i,j=l k,k'=l ^ ^ 



+ J^[^kk'GfjJco{x^-,X^-,XuJ + SijF,^f^,^J^fco{x^-,Xu^,Xu^,)] 



tanjS 



^^^xT^x^Q^j ~ ^^kQ\j) fDo{x^-,x^-,Xu^,Xy,)^ (3.41) 
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rA^Jx^fBO (x^-,Xu^^ +0{tan'^(3), 



(3.42) 



i,j=l k,k'=l '-^ ^ m=l 



K 

— -Okk' 
2 \X, 



Skk'VjiV*i 



/o^(xsJ + Go(a;^-s,,a:^-nJ 



tan 



'l3,Jx~x~U*2Ut2fD0 (a;^- , x^- , Xu^ , xc, ] 



(3.43) 



where ?7 and V are matrices which diagonahze the mass matrix of charginos, Tm (i=L, R) is the 
matrix which diagonahzes the mass matrix of the scalar up-type quarks and K is the CKM matrix. 
For the definitions of various symbols in the above equations, see the Apendex B. In eqs. ( p. 401 ) and 
( p.42| ) the tan^/5 term has been explicitly written. 
Let us give some remarks: 

(a) The first term in eq. ( p.39| ) which is propotional to sin^2a arises from the trilinear H^H^H 
couplings. At tree level, due to the more constraints ( p. 36 - ^^3^ ) than that in the 2HDM defined in 
subsection 3.1, there are only two free parameters in the Higgs sector of MSSM which we may choose 
as tan/3 and one of masses of Higgs bosons, e. g., m^o. In the large tan/? limit one has 



m 



m 



Zi 



m 



AO) 



sin2a ~ cotp. 



(3.44) 



So the first term would have no tan/3 enhencement. However, including the radiative corrections, the 
above relations are , in general, changed and the mixing angle a is determined by 



tan2a 



sin2(3{m\o + m 



2R 



12 



TT 



cos2(3{m\o - ml) + R22 - R 
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- < a < 0, 
2 - ' 



(3.45) 



where Rij are the radiative corrections to the mass matrix of the neutral Higgs bosons in the {H^, H2} 
basis and have been given in refs. pO|, [2l[| . As shown in ref.|2^, R12 can reach more than ten percents 
of R22 in the case of /i ~ ~ A;, and consequently sin 2a can be the order of one. Of course, 
there are some cases in which R12 as well as Rn is very small and of a few thousandth of R22 and 
consequently the tree level relations ( |3.44| ) are almost not changed. Therefore, we keep, in general, 
a and Higgs boson masses as free parameters because several parameters in MSSM enter the Higgs 
sector through the raditive corrections. 

(b) The first terms in eqs. ( p.40| , p.42| ) which arise from the self-energy type diagram will 
provide the tan^/? enhencement, as pointed out in [^, |^, if the mass splittings of stops are large (say, 
>100Gev). The condition is necessary because if all the squark masses are degenerate (m^-^ 



m 



t2 
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nib 








Mb, 




G p 


4.8 GeV 


1.4 GeV 


0.11 GeV 


1.78 GeV 


5.37 GeV 


0.22 GeV 


1.17 X 10-^GeV"2 




\VM 


\Vcb\ 


Br{b ceVe) 








129 


0.0385 


0.036 


0.114 


1.54 X 10-12 s 


0.232 


0.12 



Table 1: Values of the standard model parameters used in our numerical analysis. 



a 




niHO 


niAo 


rriH- 






0.08 


110 GeV 


160 GeV 


160 GeV 


180 GeV 


90 Gev 


150 GeV 



Table 2: 



m), the large contributions arising from the chargino-squark loop exactly cancel due to the GIM 
mechanism |2^. We remark that the chirality structure of the Qj(i=l, 2) operators allows a large 
tan/3 enhancement for the Wilson coefficients CQ.(i=l, 2), as happened for the magnetic moment 
operator O7, and there is no such a large tan/3 enhancement for the Wilson coefficients Cj(i=8, 9) 
due to the different chirality structure of the Oi(i=8, 9) operators. 

(c)The last term in eq. ( 3.43 ) comes from the chargino-chargino box diagrams and are propor- 
tional to tan^/?. It was missed in the literature. For l=/i, r and large tan/5, it is numerically the same 
size as the other contributions so that it should not be omitted. 



4 Numerical results 

Below we assume no CP violating phases from 2HDM and SUSY. As said in section II, there are 6 
free parameters in the 2IIDM which are tanf3, a, rriho, mno, rriAo, mH±- In MSSM, in addition to the 
above 6 parameters, 7 extra free parameters, mg^, m^^, m^^. At, M2, /i, and m^;^ are needed in order 
to calculate the Wilson coefficients. In Table. 1 we list all SM inputs for our numerical analysis. 

Numerical results are given in Figs. 3-9. Figs. 3-5 are devoted to the decay Bg — > in 
2HDM. In the numerical calculations in 2HDM the constraint on m/^± from b —* 57 1^ has been 
imposed. We present in Fig. 3 the branching ratio (Br) for Bs fJ'^fJ'^ as the function of mfjo, 
the mass of the heavier CP even neutral Higgs boson, for fixed values of the other parameters. The 
figure shows that the Br increases when mno increases except for the mixing angle a=0. The reason 
is that in the large tan/5 case the trilinear H^H^H {B.=H^, h^) couplings are proportional to m'j^o 
and/or m\o and the couplings vanish when a is equal to zero. One can see from Fig. 4 that the Br, 
as the function of a, behaves like sin22a when mfjo is large enough ( say, 500 GeV), which implies 
the contributions from the trlinear couplings of NHBs dominate (see eqs. (2.13), (3.30) and (3.31)). 
Fig. 5 shows the tan/3 dependence of Br and one can see from it that the contributions coming from 
Cq^ s can dominate when tan/? is large enough (say, larger than 80). In most of the large tan/? region 
in the parameter space Br is of order 10~^, an order of magnitude larger than that in SM. 

The numerical results in MSSM are presented in Figs. 6-9. We present the correlation between 
C7 and for I = fi and I = t respectively in Fig. 6 and 7 where the absolute values of Cj are 
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taken from the data of B ^ Xs'~f[2^ with the 2a errors imposed. We set m^^, m^-^, m^^, rn^-, 
and tanf] as random free parameters. They vary in the range 180-300 GeV, 250-450 GeV, 200-400 
Gev, 160-360 Gev and 25-50 respectively. Other parameters are fixed as given in Table 2. We get 
about 3000 permitted points among 25000 points. The contributions to Wilson coefficients due to 
superparticles in a loop ( SUSY contributions) come mainly through the ulXid vertex, which is 
determined by the mixing between Higgsinos and Winos and the mixing between stops. The vertex 
appears in Feynman diagrams which describe the processes b ^ s'-f and b sl~^l~ so that there exists 
a correlation between C7 and Cq-. In some large tan/? regions of the parameter space in MSSM, 
SUSY contributions interfere destructively with the SM contributions and SUSY contributions can 
be so large that they can overwhelm those from the SM and the Higgs sector so that the sign of C7 
is changed compared to that in SM, as shown in Fig. 6 and 7. In the regions Cq.s are proportional 
to tan^/3 and consequently can compete with Cg for l=/i and be much larger than Cg for l=r. We 
also calculate the correlation between Cg^ and Cq^ in the regions and it follows that Cq-^ ~ —Cq^. 

The contribution to Cg which is proportional to tan^jS coming from chargino-chargino box dia- 
grams is numerically the same order as other contributions from chargino-chargino box or chargino-up 
type squark penguin diagrams. As a whole SUSY contributions to Cg give about 10% corrections 
to the SM value. Taking Cq^ in the allowed range in Figs. 6 and 7 and Cq^ ~ — Cq^, we draw the 
branching ratio of Bg — > as function of Cq^ in Figs. 8 and 9, given Cg being the SM value with 
10% variations. Figs. 8 is for l=/i and Fig. 9 for l=r. From Figs. 8,9, we can see that Br(5s — >■ /"'"/") 
is more sensitive to Cq_ s (which represent the contributions from NHBs ) than to Cg in the large 
tan/3 (larger than about 30) regions because in the regions |Cq.'|s are much larger than |^^^Cg|. The 
numerical results tell us that it is possible to saturate the experimental bound (1.1) for Bg /i^/i^ 
in some regions of the parameter space in MSSM. In other words, the experimental bound could 
impose a constraint on the parameter space of SUSY which we shall analyze elsewhere. Because 
Cq^ s are proportional to the lepton mass, Br for l=r can reach order of 10~^ in the regions in which 
Br(i?s /W^A''") saturates the experimental bound. 



5 Conclusions 

In summary we have analyzed the decays B^ — > in the model II 2HDM and SUSY with large 
tan/3. Although these decays have been studied in these models before and reanalyzed recently, it 
seems that no complete analysis exists so far. We have calculated all leading terms in the large tan/3 
limit. We found that in addition to the Higgs boson- W boson box diagram, the chargino-chargino 
box diagram gives also a contribution proportional to tan^/3, the former to Cq- (i=l,2) and the latter 
to Ci (i=8,9). The contributions from NHBs always increase the branching ratios in the large tan(3 
case so that the branching ratios in the 2HDM and in SUSY are larger than those in SM. We have 
numerically computed the branching ratios for l=/i and r. In the 2HDM the branching ratio for 
l=/i is about 10~®, an order of magnitude larger than that in SM, if tan/? = 50 or so and the other 
parameters are in the reasonable range. We have shown the dependence of the branching ratio with 
respect to the mixing angle a and neutral Higgs boson masses. The branching ratio increases when 
the splitting of the masses of the two CP even neutral Higgs bosons increases except for the case of 
the mixing angle a=0. In MSSM the branching ratio for l=/i can saturate the experimental bound 
in some regions of the parameter space where Cq.s (i=l,2) behave as tan^/3. In the other regions 
where Cg.s (i=l,2) behave as ia.v?(3 the branching ratio is about the order 10~^. The branching 
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ratio for l=r reaches 10 ^ in the regions of the parameter space in which Bt{Bs l^l ) saturates 



the experimental bound. In the near future when very high statistics can be reached p|, ^ the 



measurements of the decays Bg l^l (l=/i, r) can provide a large potential to find or exclude the 
large tan/3 parts of the parameter space in 2HDM and/or SUSY. 
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Appendix 

A Wilson coefficients in a model II 2HDM 

Wilson coefficients are extracted from the transition amplitudes by integrating out heavy particles. 
In our convention, the effective Hamiltonian is related with the amplitude by 

iM{b -> sl+r) = -i< l^rs\Heff\b > (5.46) 

By computating the self-energy type, Higgs-penguin and box diagrams, Cg^ and Cq^ with the su- 
perscript denoting the type of a diagram are extracted out, as given below 

Cqi = facrhnixH- - l)xtf2ixt,yt) , (5.47) 

Cq2 = -facTAixH- - l)xtf2ixt,yt) , (5.48) 

P _ sinHa {mlo - mlo)\ ^ 

- fac [-1 + (xh- - l)rhH] Xtf2iXt, Vt) , (5.49) 

Cq2 = fcic [-1 + {xh- - IYa] Xtf2ixt, Vt) , (5.50) 

= -facB+{xH+,Xt) , (5.51) 

= facB+{xH+,Xt) , (5.52) 



where 



^ .-o.-t -.n lyrz I sin^ a cos^q;\ ,_ 
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l-y + ylny 
•^i(^) = (^_i)2 ' (5-55) 

f2{x,y) = - ^-5;^^ — - + - ^ -, with z = x/y, (5.56) 

[z — x)[x — 1) [z — l)[x — 1) 

B+{x, y) = , 5.57) 

X — y y — 1 X — 1 

B Wilson coefficients in MSSM 
B.l Feynman rules and conventions 

In this subsection we present our convention. In order to avoid the trouble in dealing with charge 
conjugate operation, we choose x" ^ the particle. The interactions of duXi HxXi Huu, and Hdd 
can be expressed as: 

^dux — '^^IXi i^ijk Pl + 'Hijk PR\dj + h.c. , 
92 



^Huu = -^j^ fijh ulUjHl^ + h.c. , 
^Hdd — Qh ddH^ , 



where 



r,,, - -,/2V,,{Tu,K),, + V,,^—^ , 

^ _.r iTuLKmd)kj 
- Mwcos(3 ' 

g'lj^ = V2 {{Qijsma - Sijcosa), -{Qijcosa + Sijstna), -t{QijSin(3 + Sijcosf])},^ 
g^jf^ — ^(QljSina — SjjCosa), —(QljCosa + S\jSina),i{Q\jSin(5 + S\jCos(5)^^ , 

rud fSina cosa . 

9h = { ^, -,itanp-f5}h , 

V2Mw cos(3 cosfJ 

fijh — TiiHl^iTiij , Qij = ---^Ui2Vji , Sij = ---^UiiVj2 . 
In (7.58) and (7.59) K is the CKM matrix, and m„ and are defined as 

■ifiu = diag{mu, rric, , rhd = diag{md, rris, nib} ■ 
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can be expressed as 



* \_sinp Mw cos-^9w 



V2Mm 



V sinp " smp 



cosa rh^ _ sin{a+p) r-j 



V2Mw 



sm/3 Mw ' cos^tlu] 



sinfd ^ 





x/2 



(- 



Inl + SHt«^t 

\/2Mi4/ V sin(3r^ sinfi 

cos-^9w 



sin/3 M14/ 



—im. 



where 



-D„L = (o - euSin^Ow) , Dur = eusiv?6^ 



The 6x6 mass matrix of u-type squark is given as 



( ( 9 ~ 9 
{m% + m 



Q 



+ cos2/3DuLml 



'-ix*ctg(5l + A) rhu 



m 



m 

2 I ™2 



(-ixctg[3l + 



(5.69) 

(5.70) 
(5.71) 

(5.72) 
(5.73) 



where each block is a 3 x 3 matrix. A is defined by Y"^ = Ak^, while is the trilinear coupling 
matrix of up-type squarks, and is the Yukawa coupling of up-type quarks. The 6x67 matrix is 
defined as 



ki'~i' 

UR^R- 



The convention of chargino masses is given as 



X 



Ms V2Mwsinl3 ) 



where 



, V2Mwcos/3 n J ' 

U*XV^ = diag{m ±,m +} 



With the above conventions, it is straightforward to extract Feynman rules. 



(5.74) 

(5.75) 

(5.76) 

(5.77) 
(5.78) 

(5.79) 



B.2 Cg, Cq, and Cq, in SUSY 

Wilson coefficients are extracted from the transition amplitudes by integrating out heavy particles. 
Cg is given as 



C9 — C9,z + Ci 



9,z 



^9,z + ^9,z + ^9,z 



(5.80) 

(5.81) 
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c, 



9,B 



Asin^Ow 



1 



(5.82) 
(5.83) 

(5.84) 



11^® r 3 



— 



-[G{xtw,0)-G{Q,Q)] 



a 



w 



1 1 



(5.85) 
(5.86) 



2 6 



9 ^^jbk^^ 
W ^ ij=i k=l 



isk 



^ ton^/J^x^- x^- U*2Ui2fD0 {x^- , x^- , Xu, , x,^ ) 



(5.87) 



By computating the self energy type, Higgs-penguin and box diagrams, Cq^ and Cq^ with the first 
superscipt denoting the type of a diagram and the second superscript a Higgs boson or a superparticle 
in the loop of the diagram are extracted out, as given below 



_ ss 

Ql(2) '^Ql(2) '^Ql(2) 



^Ql(2) 



riPH I r<PS 

^Ql{2) ^ ^Qn2) ' 



^Q2 



riPH 



-facThHixH- - l)xtfc°{xH-A,Xt) , 
facTAixH- - l)Xtfco{xH-,l,Xt) , 

sin'^2a (m\o — tn^oY 



(5.88) 
(5.89) 
(5.90) 
(5.91) 
(5.92) 



-VthiVt) 



2 2m]jom\o 

+fac [-1 + (xh- - l)rhH] Xtfc^ixH-, 1, Xt) , 
-fac [-1 + {xh- - 1)^^] Xtfco{xH-A, Xt) , 
Cgf = fac [fco{xH-,l, Xt) - fco{xH-,l, 0)] , 
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^Q2 



(5.93) 
(5.94) 
(5.95) 



C|f = -fa.[fc'>{x„-,l,Xt) - fco{xH-, 1,0)] , 



2 6 
" * i=l fc=l 



(5.96) 
(5.97) 



2 6 



(5.98) 



— /a 



w 



2 6 



mbtan(3X- 



isfe 



* i,j=l k,k'=l 



(5.99) 



2 6 



'^Oi - /aczrT735^ 2^ riihk'^ isk 



rribtanPX- 



* i,i=l k,k'=l 



^kk'GfjJco {x^- , x^- , Xu^] 



riBS 



'fa 



2 6 



(5.100) 



w 



mi,tan(3\t 



) (5.101) 

,j=lk=l \v ' ^ / '3 



fac ^ E E (^^x-^x-<5u - 2;«fc<5l j /^°(^xr ' ^xj ' (,5.102) 



where 



1 mfem^ 2 2 



/ sin^ a cos^ a' 
+ 



TA = , Xh- = 



2 



m 



AO 



) — 



m. 



1 



(5.103) 
(5.104) 



^ijk ™2 



+ - ^ 



sma[ . /a^^-a^^-^fto + ^^uQiju^) 



m 



2^^cosa{^^I~I~gl.jjo+Xu^g..Ho) , (5.105) 



^ j^2^ 
^Ufe = ^^^^\/^X-^X-3ijA0 + ^Uk9ijAo) ' 

^kk' = — ^sinajkk'h° ^cosaJkk'RO , i'kk' = t'^rhk'AO ■ 



m 



m 



HO 



m 



(5.106) 
(5.107) 



AO 



The related loop-integral functions are defined as 



F.(x) = 



2(x - 1)3 



x^ -4x + 3 + 2logx 



(5.108) 
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FAx) 



Fo(x) 



1 



2(x- 1)3 

1 

2(^-1)2 



Fo{x,y) 
Goix,y) 

G'{x,y,z) 
fBoix,y) 
fco{x,y,z) 
fDo{x,y,z,w) 



x-y 
1 

x-y 
1 

x-y 
1 

x-y 
1 

x-y 
1 

y- z 
1 

z — w 



X 



X — 1 
.2 



X — 1 — 2x logx 

— 7x + 6 + (3x + 2) logx 
logx — {x ^ y) 



OC^ 3 

logx X — (x ^ y) 

X — 1 2 



x^ 1 

logx — {x ^ y) 



X — 1 X — 1 

[Goix,z) - Go{y,z)] , 

{xlogx - ylogy) , 

ifBoix,y) - fBo{x,z)) , 

{fco{x,y,z) - fco{x,y,w)) . 



5.109) 
5.110) 
5.111) 
5.112) 

5.113) 
5.114) 
5.115) 
5.116) 
5.117) 



In the above expressions for Wilson coefficients, the possible mixing among sneutrinos has been 
neglected. 
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Figure Captions 

Fig. 1 Feynman diagrams which give the leading contributions to Cg^ and Cq^ in the 2HDM 
with large tan/3. 

Fig. 2 Dominant Feynman diagrams in which the virtual superparticles are exchanged in the loop 
in MSSM with large tan/?. 

Fig. 3 The branching ratio of Bg ^^^^ as functions of Mfjo in the 2HDM. Curves labelled by 
1, 2, 3 corresponds to a = 0,7r/8, 7r/4 respectively. Other parameters are chosen to be Mf^o = 120 
GeV, Mh± = 250 GeV and tanP = 60. 

Fig. 4 The branching ratio as functions of a in the 2HDM. Four curves labelled by 1, 2, 3, 4 
corespond to Mho = 220, 320, 420, 520 GeV respectively. Other parameters are the same in Fig. 3 

Fig. 5 The branching ratio as functions of tanjS with a = 7r/8 in the 2HDM. Four curves are 
classified as in Fig. 4. Other parameters are chosen as in Fig. 3. 

Fig. 6: The correlation between and Cq^ for l=fi in MSSM with mj^, m^-^, mg^, m^- and tan(3 
as free parameters. 
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Fig. 7: The correlation between C7 and Cq^ for \—t in MSSM with mf^, mf^, fric^^ m^- and ton/3 
as free parameters. 

Fig. 8: Br{Bs f^^fJ^~) as a function of Cqi in MSSM. Cqi a; — Cq2 has been assumed. Three 
hues correspond to Cg = —4.2, —4.6, —5.0 respectively. 

Fig. 9: Br{Bs r+r~) as a function of Cqi in MSSM. Cqi — Cq2 has been assumed. Three 
lines correspond to Cg = —4.2, —4.6, —5.0 respectively. 
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Figure 1: 
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